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DIFFUSION LIMIT FOR MANY PARTICLES IN A 
PERIODIC STOCHASTIC ACCELERATION FIELD 



By Yves Elskens and Etienne Pardoux 

'do' 

I CNRS-universite de Provence 



The one-dimensional motion of any number A/" of particles in the 
field of many independent waves (with strong spatial correlation) is 
^ ' formulated as a second-order system of stochastic differential equa- 

O I tions, driven by two Wiener processes. In the limit of vanishing parti- 

. cle mass m ^ 0, or equivalently of large noise intensity, we show that 

the momenta of all TV particles converge weakly to A/" independent 
Brownian motions, and this convergence holds even if the noise is 
periodic. This justifies the usual application of the diffusion equation 
' to a family of particles in a unique stochastic force field. The proof 

, rests on the ergodic properties of the relative velocity of two particles 

T ' in the scaling limit. 

1. Introduction. Th e motion of a particle in the field of many waves 
DoGuOi ndamental process in classical physics, the 



understanding of which is a prerequisite to the analysis of many plasma 
and fiuid phenomena. In one space dimension, it can be described by the 
' Hamiltonian model 

g. p2 M 

(1.1) H = n h V Am COs{kmq - UJmt - ^m) 

771=1 



where the particle with mass m has position q and momentum p, while the 
force field derives from a potential with time Fourier components Amc"^™. 
The wave field comprises M waves, with a smooth dispersion relation asso- 
X ■ elating a wavenumber km , sl pulsation and a phase velocity Vm — i^m/km 

■ to each wave ~ usually determined by fixed properties of the environment, 

such as the geometry of the domain where waves propagate (then wavenum- 
bers km and pulsations 

ujm are discrete) . The complex amplitudes AmS^'^"^ 
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are more easily tuned by the experimenter or affected by simple changes in 
the environment. 

The dynamical systems approach to this problem discusses the particle 
motion after prescribing a single choice for each wave complex amplitude. 
As it would be quite exceptional to control all wa ves (thou gh this is e.g. 
the assumption underlying the standard map, see BeEs98bl | for a discus- 
sion), physicists often turn to a probabilistic description of the dynamics, 
considering an "ensemble" of realizations {Am,^m)- Various arguments are 
then invoked to reduce the particle evolution equations to a stochastic dif- 
ferential equation, often dr iven b y a "white noise" . This results in somewhat 
tractable models (see e.g. Kr02| about the validity of such derivations). 

In this paper we focus on two issues. First, a random field characterized by 
(Amjifrn) for a given dispersion relation with discrete frequency spectrum 
may be periodic in time : may the force on the particle be considered as 
independent over several time periods in a genuine limit ? Second, may one 
consider several particles subject to the same wavefield as independent in 
a genuine limit ? The latter issue underlies the frequent application of the 
Fokker-Planck equation to the evolution of a family of particles in a single 
turbulent wavefield - though a priori one can only grant that the diffusion 
equation describes the evolution of the distribution of a single particle for 
an ensemble of wavefield samples. 

From a more general perspective, this work also r elates to the issue of 
"propagation of chaos" in statistical physics KaSd . Ka59l |. an aspect of 
Hubert's 6th problem : how does chaotic dynamics enable a system, in which 
initial data are independent ("random") but the evolution may generate cor- 
relations, to behave as if the evolution regenerated independence ("random- 
ness") or destroyed correlations ? Here, how do two Wiener processes, fully 
describing a prescribed "turbulent" environment, generate independent 
Brownian motions for particles ? 

A further motivation for the present work is that physics literature most 
often focuses on the evolution of particle distribution functions, e.g. by show- 
ing that they obey a Fokker-Planck equation, and on instantaneous observ- 
ables such as p{ti) for given ti (pointwise in time). However, the notion of 
a diffusion process implies rather a measure on the set of trajectories, viz. 
functions p{-) (globally in time). Here we shall show how our model im- 
plies that an arbitrary number of trajectories in a single realization of the 
dynamics do, jointly, admit the Wiener measure description. 

In sections [2] and [3] we motivate the mathematical model more precisely 
and state our main results, which are proved in the subsequent sections. 
The crucial Theorem 13.21 is an ergodic theorem for a rescaled process, im- 
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plying that a process Vj, describing the relative velocity of one particle with 
respect to another one (or to its own earlier motion), converges weakly to 
a Brownian motion. The difficulty in proving the ergodic theorem is that 
the invariant measure of the diffusion process is infinite (it is the Lebesgue 
measure dxdy), and we must estimate a continuous additive functional gen- 
erated by a function (namely f{x,y) = sin^x) which is not integrable with 
respect to this measure (it is only locally integrable). This weak convergence 
then implies the final many-particle result (Theorem 13. ip by a straightfor- 
ward application of the Levy characterization of Brownian motion. Section 
[8] outlines implications and possible extensions to this work. 



2. Physical background. Because the waves have different frequen- 
cies and velocities, it is generally unrealistic to assume their phases to be 
correlated. Their intensities are more easily observed, but both in nature 
and in the laboratory the accumulation of statistical data on waves often 
involves only their average power spectra, not the detailed intensity data 
for each measurement run. We assume here that these complex amplitudes 
are random data, and investigate the statistics of the particle motion in 
the resulting time-dependent random field. This dynamics is a "stochastic 
acceleration prob l em" for a "passive p article" in weak plasma turbulence 
DiKrSd boGrSl iMaElsd . ISted . IVEQ?! . and its understanding is a prereq- 



uisite to a proper anal ysis of the case w here the particle motion feeds back 
on the wave evolution |DoCa97 . ElEs03 |. 

The Hamiltonian (jl.ip generates equations of motion 



(2.1) 
(2.2) 



p/m, 



m 



An important observation Chi79l . Es85 | on the motion of a particle in 
the field ()2.2p is locality in velocity : the evolution of the particle when it 
has velocity q = v depends only weakly on the waves with a Doppler-shifted 
frequency cOm — k^v much larger than their trapping oscillation frequency 
km^/ArnJm. In particular, for a two- wave system the resonance overlap pa- 
rameter 



Sl,2 



\iV2/k2 - toi/ki 



becomes unity when there exists a velocity u 
2^JAi/m (with fci > 



= u:2/k2-2^/M/m = uJi/ki + 
For many waves with overlap 



0, k2 > ijJ2 > OJl) 

parameters s ^ 1, the relevant phase velocity range for waves influencing 
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the parti cle is a "resonance box", with a width scaling as {A/m)^/^ . |BeEs97 . 
BeEs98al | 

A good approximation to typical wave dispersion relations in the strong 
overlap limit, after a Galileo change of reference frame (see e.g. sec. 6.7 in 
Hisoi), is 



(2.3) 



km = ko, ujm = 27r(m - M/2)/T, 



for some /cq, T. Then, in the limit A4 — > oo, the equations of m otion vield 
for Ame^^"^ = Aq (with real ^o) the well-known standard map BeEs98b |. 
The case where phases ipm are independent random variables uniformly dis- 
tributed on the circle [0, 27r], while Am = Ap is given, was investig ated 
notably by Gary, Escande, Verga and Benisti |BeEs97l . lGEV9nl . IeIEsO.^I I and 
occurs in the context of the random phase approximation. 

To the extent that the phases and amplitudes of the waves are independent 
random variables, the physicist usually views the force ()2.2p as a mollification 
of a white noise, with amplitude a = ^/K{k':^Af^) = \km\\/^^A^, where the 
relevant mode m is the one nearest to the current particle velocity (the 
mathematical expectation IE is called ensemble average with respect to wave 
amplitudes and phases) 



This is the core of quasilinear theory [DrPi62l . |Pe94l . iRoFiGll . |VVS62| |. 
With some care, one interprets (I2.ip - (l2.2p as a stochastic diff'erential equa- 
tion ; this applies in the case m G Z for dispersion relation (j2.3p with Gaus- 



sian independent complex amplitudes such that KAm = k^ ^u^. The particle 
velocity then has a Brownian evolution, so that for t,t' £ [0,T] 



(2.4) 



^{Pt - Po){Pt' - Po) = -Dql min(t, t') 



with the quasilinear diffusion coefficient -Dql = However, the par- 

ticle evolution for t > T may show a strong correlat ion t o its m otion for 
< t < T because the waves are periodic in time ElOTI . |E108| |. and the 
dispersion relation (j2.3p may generate a strong spatial correlation between 
the motions of two particles because all waves acting on a particle at any 
time have the same wavelengthH 



^Phases do not appear in a (nor in s) because the relative phase of two waves ifm + 
ujmt — fn— ^nt varies uniformly over time (hence iprn — 'fin can be absorbed in the choice 
of the time origin) . 

^ There is a large body of literature on the case of incoherent waves with no dispersion 
relation. Then the sum becomes a double sum and one varies wavenumbers kn 

independently from pulsations ujrn- This space-time stochastic environment is more noisy 
than our model and may also be considered to motivate a quasilinear approximation. 
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Set ko = 1 and T = 27r by the choice of space and time units. The 
particle phase space is T x M, where T is the circle modulo 2tt, and the 
particle equation of motion reads, with initial data {q{0),q{0)) = {qo,qo), 

q = ^ Am sm{q - mt + fm) 

(2.5) m=-fi 

= ^ Am cos(mt — ipm) sin q — m^^ ^ Am sm{mt — ipm) cos q 

m m 

in a galilean frame moving at a velocity inside the spectrum of wave phase 
velocities {fi + fi' = Ad). We assume 3> 1, /i' » 1. For Am = Aq given and 
independent random phases (uniform on the c ircle) , in the limit An /m oo, 



with {m/Ao)'^^^M > 10, Benisti and Escande |BeEs97l . lBeEs98al | have shown 
that the particle momentum p = mq follows essentially a Brownian motion, 
with diffusion constant given by the quasilinear estimate 

Dql = ttAI 

as long as the motion does not approach the boundaries of the wave velocity 
spectrum. Then the particle momenta p for an ensemble of independent re- 
alizations of the system will be described by a distribution function verifying 
the Fokker-Planck equation 

(2.6) dtf = ^dlf 



even for t > T. Numerical simulations E108l | show similar behaviour for i.i.d. 



complex Gaussian random variables Am^^'^'"- The present work establishes 
a rigorous version of this result in the frame of stochastic processes. 

3. Main results. We first let min(/i, /x') ^ oo in the model, taking 
Am cos ifm, Am sin ifm as i.i.d. Gaussian random variables with zero expecta- 
tion and E,(Am cos (fm)'^ = ^{^m simpm)'^ = 1/2 for all m G Z. In particular 
this implies that phases <fm are i.i.d. uniformly on T. To follow earlier prac- 
tice, we now set m = 1/A. Then formally ()2.5|) becomes the Stratonovich 
stochastic differential equation for < i < 27r 

(3.1) dQt = APtdt, 

(3.2) dPt = sm{Qt) o dCt + cos{Qt) o dSt , 

with initial data Qo = qo, Po = Po = ^Qo = Qo/A ; here, from KaSSi ]. 



vr ^/^(C, 5) is a standard 2-dimensional Brownian motion. In other words, 
C and S are martingales, with Cq = Sq = and 

(3.3) {C)t = {S)t = TTt , {C,S)t = 0. 
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Note that the Stratonovich and Ito integrals define the same solutions for 
this system, and that the vector fields {sin q)dp and (cos q)dp commute. 

For < t < 27r, it is clear that P is a Brownian motion for any value 
oi A > (see the first lines of Sec. [7]). However, the model (|2.5p defines a 
dynamical system for all times < t < oo, and one may wonder how the 
initial stochastic behaviour over [0, 2tt] extends for longer times. Formally, 
one solves then (|3.ip - (|3.2p with the periodized field, i.e. with the continuous 
processes defined by 

(3.4) dCt+2k7T = dCt , dSt+2k7T = dSt 

for A: G Z. In other words, Ct — ■^C2tt and St — ■^S2tt are independent 
Brownian bridges repeated periodically for t G R, while C2tt and S2-K are 
independent Gaussian random variables with expectation and variance 

For this extended process, the wave field acting on the particle for t 
[0, 2tt] is not stochastically independent from the wave field acting during 
[0, 27r]. Therefore one does not expect the particle momentum to proceed as 
a Brownian motion for all times, and indeed for A small enough the velocity 
q may remain bounded in a narrow interval for all times. This is easily 
seen numerically and can be attributed to the existence of Kolmogorov- 
Arnol'd-Moser invariant tori in the 3-dimensional extended phase space with 
coordinates {t,q,q). 



On t he other hand, for large A, the dynamics viewpoint |BeEs97l . lBeEs98al . 



suggests that the nonlinearity in the equations of motion (due to 



trigonometric functions of Q) may enable a decorrelation of the force over 
the period T = 27r, so that the long-time evolution of the velocity would 
also be close to Brownian. This is what we shall show. 

An intimately related issue is the relative motion of several particles, re- 
leased in the same realization of the wave field. Even though each particle 
velocity diffuses for t £ [0,27r], their motions are not independent. We shall 
also show that for large A the motions of any finite family of particles re- 
leased at initial data {Q^q\ Pq'^^), I < v < M, approaches a family of J\f 
independent processes. This can also be expected from the consideration of 
the top Lyapunov exponent of the dynamics (|3.ip - (|3.2p in the limit A oo. 



Theorem 3.1. For any M > 0, the momentum processes P^"^ defined 

by 

(3.5) dQS^) = AP^Ut, Qt^=q^^\ 

(3.6) dP/") = (sinQS"))dCt + (cosQS'^))d5t, P^""^ = vt\ 
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withM different initial data {Qq^Kpo'^) ^ TxR, 1 < v < J\f , converge as 

oo to M independent Wiener processes with variance nt, and convergence is 

in law in C{M.+,M.^). 

The key argument in the proof is the following weak convergence theorem, 
where we write now n = n^^'^A. Consider the two-dimensional diffusion 
process indexed by n > 1, solution of the SDE on 

( dur 



(3.7) 



nVr,Uo = u, 
dV^i" = sm{U{')dWt , Vq = v 



where {u,v) {(/cvr, 0), k £ 7j} and PV^ is a standard Brownian motion. We 
prove 

Theorem 3.2. As n ^ oo, 

where {Bt, t > 0} is a standard one-dimensional Brownian motion, and the 
convergence is in law in C(M+,M). 

4. A change of time scale. Note that for any n > 1, the law of 

{(f7", Vf^), t > 0}, the solution of (|3.7p . is characterized by the statement 

-^ = nVr, Uo = n, 

is a martingale, ^^^^ = sm'^iUD , Vq = v ■ 
Now define (hke Benisti and Escande BeEs97 . BeEs98al |) 
(4-1) ^* = f^:-/3,, Y, = n'/'v:.y,^. 

We first note that Xq = u, Yq = n^^^v, y is a martingale, and 

f^ = --^/^^(--^/^t)=nV3y:„.,3, = r., 

Using a well-known martingale representation theorem, we can pretend that 
there exists a standard Brownian motion {Bt, t > 0} such that 

(4.2) ^^ = ^*'^° = "' 

dy* = sin{Xt)dBt , Yo = n^l^v . 
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Note that the process {{Xt,Yt), t > 0} stih depends upon n, but only 
through the value of Yq. 

On the other hand, V^^ = n~^/^y^2/3j. Hence 

Vt'' = v + / sm{Xs)dBs , 

Jo 

in other words V"' is a martingale such that Vq* = v and 



Jo 



Here we recall the fact that the process X depends upon n (through the 
initial condition of 1"), unless v = 0. Consequently 

(4.3) lim =tx lim — - / sm'^(X^)ds 

and in order to prove Theorem 13.21 it suffices to show that the above limit 
is t/2. 

5. Qualitative properties of the solution of ()4.2p . We now consider 
the two-dimensional diffusion process 

r dXt 



(5.1) 



Yt, Xo = x, 



dt 

^ dYt = sm{Xt)dBt , Yo = y, 

with values in the state-space E = [0,27r) x M\{(0, 0), (vr, 0)}, where 27r 
is identified with 0. We first prove that the process {{Xt,Yt), t > 0} is a 
conservative S-valued diffusion. Indeed, 

Proposition 5.1. Whenever the initial condition {x,y) belongs to E, 
inf{t > 0, {Xt,Yt) G {(0,0), (^,0)}} = +oo a.s. 

Proof: We define the stopping time 

r = inf{t, {Xt,Yt) = (0,0)} . 

Let Rt = Xf + Y^^, Zt = log Rt, t > 0. A priori, Zt takes its values in 
[— cxD,+oo). Ito calculus on the interval [0, r) yields 

dX^ = 2XtYtdt, 

dY^^ = 2YtsmiXt)dBt + sm^iXt)dt, 
di?t d{R)t 

2YtXt + sin^ (Xt) , Y.^ sin^ (Xt) , Yt sin(Xj ) , „ 
= ^—^dt - 2-^ j—^dt + 2 ^ dBt . 

Rt Rt Rt 
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Now clearly | sin(x)| < |2;|, sin^(x) < x^, and it follows from the above and 
standard inequalities that on the time interval [0,r), 

Zt>Zo-2t+ f ifsdBs , 
Jo 

where l^^^l < 1. Hence the process {Zt, t > 0} is bounded from below on any 
finite time interval, which implies that r = +oo a.s., since r = inf{t, = 
— oo}. A similar argument shows that r' = +cxd a.s., where 

T' = inf{t, iXt,Yt)e{{0,0),{7T,0)}}. 

□ 



We next prove (here and below Be stands for the u-algebra of Borel 
subsets of E) the 

Proposition 5.2. The transition probabilities 

{p{{x, y); t, A) := W,J{Xt, Yt) G A), (x, y)eE, t>0, Ae Be} 

have smooth densities p{{x,y);t, (x' ,y')) with respect to Lebesgue's measure 
dx'dy' on E. 

Proof: Consider the Lie algebra of vector fields on E generated by Xi = 
sin(x)^, X2 = [Xo,Xi] and = [Xq, [Xq, Xi]], where = y^. This 
Lie algebra has rank 2 at each point of E. The result is now a st andar d 
consequence of the well-known Malliavin calculus, see e.g. Nualart Nu06l |. 
□ 



Proposition 5.3. The E~valued diffusion process {{Xt,Yt), t > 0} is 
topologically irreducible, in the sense that for all (x, y) £ E, t > 0, A £ Be 
with non empty interior, 

V^,,y{{Xt,Yt) e A) > . 
Proof: From Stroock-Varadhan's support theorem, see e.g. Ikeda-Watanabe 



IkWa89f | ■ the support of the law of {Xt, Yt) starting from {Xq, Yq) = {x, y) is 
the closure of the set of points which the following controlled ordinary differ- 
ential equation can reach at time t by varying the control {u{s), < s < t} 
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in the class of piecewise continuous functions : 



(5.2) 



' dx 

— {s)=y{s), x{0) = x; 
as 

dy 

— (s) = sin(x(s))u(s) , y{0)=y. 
as 



It is not hard to show that the set of accessible points at time t > by the 
solution of (15. 2p is dense in E. The result now follows from the fact that the 
transition probability is absolutely continuous with respect to Lebesgue's 
measure, see Proposition 15. 2[ □ 



We next prove the 
Lemma 5.4. 

F{\Yt\ ^ oo, as t ^ oo) = 0. 
Proof: The Lemma follows readily from the fact that 

Yt = W (^j\m''{Xs)dsy 

where {W{t), t > 0} is a scalar Brownian motion. Then either /q sin^(Xs)ds 
is bounded and It is finite, or the integral diverges and Yj is finite anyway 
because W is recurrent. □ 



Hence the topologically irreducible £^-valued Feller process {{Xt, Yt), t > 
0} is recurrent. Its unique (up to a multiplicative constant) invariant measure 
is the Lebesgue measure on E, so that in particular the process is null- 
recurrent. It then follows from (ii) in Theorem 20.21 from Kallenberg Ka02l | 



Lemma 5.5. For all M > 0, as t ^ co, 



1 

~t Jq 



l{|y,|<M}ds a.s. 



6. A path decomposition of the process {(Xt,lt), t > 0}. We 

first define two sequences of stopping times. Let Tq = and 

for £ odd, = inf{t > T^^i, \Yt\ > M + 1} , 
for i even, = inf{t > T^.i, \Yt\ < M] . 
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Let now tq = Ti. We next define recursively {rfc, k > 1} as follows. Given 
Tk-ij we first define 

Lk = sup{e > 0, rfc_i > T2£+i} . 

Now let 

[T2Lk+3 if Tk-l > T2Lk+2 ■ 

We now define 

Tfc = inf{t > r]k, \Xt -Xr,J= 27r} A inf{i > r]k, \Yt -YrjJ>l}. 
It follows from the above definitions that 

a statement which will be refined in the proof of Proposition 16.31 Define 



ft °° fTkAt ft 

/ Ijiy i>M+ijsin^(-'^s)ds < V / sin^(Xs)ds < / sin^(Xs)ds, 

Jo rr^ JvkM Jo 



We first prove the 
Lemma 6.1. 



K' = {k>i, -y,j<i}, 
K' = {k>i, |y.,-y,j = i}, 

Kt ={k> 1, rik<t}, 



in as M ^ oo, uniformly in t > 0. 

Proof: We shall use repeatedly the fact that since | > M > 2, jy^^. | — 1 > 
|yr)ft|/2. We have that (see the Appendix below), since — rjk < 47r/|y^j.|, 



Gi^'l^,,,) < p sup |yt-y^j>i 

\Vk<t<Tk 



■^Vk 



< 2exp(-|y^J/(8^)). 
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Consequently, using again the inequality Tk — r]k < 47r/|y^j. |, we deduce that 

Stt 



E 



{Tk- Vk)'i-{keK^}\^Vk 



< 



< 



\Yr 



Vk I 

Stt 



exp(-|y,J/(87r)) 
exp(-M/(87r)) . 



'nk\ 



On the other hand, whenever k e 

rfc-%>27r/(|i;,J + l)>7r/|y^J 
Now, provided t > Att/M, 



2t>t + 



M 



> E 



{Tk - Ilk) 

kGK° 



> ttE 



> -E 

- 2 



J2 MkeKoyjTri 

k€Kt ^ 



k&Kt I 



Since 



¥{k G i^°|^r,J = 1 - P(fc G K^\Tr,J 

> l-2exp(-M/(87r)) 
>l/2, 



provided M is large enough. Finally 
1 



-E 



{Tk - Vk) 

keKl 



< 32exp(-M/(87r)) 

= 32exp(-M/(87r)) 
^0, 



^[EkeKM-'] 



as M — > oo, uniformly in t. 



□ 
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Now, for any k G K^, 



Tk 



S\T?{Xs)ds 



Vk 



2tt 



(6.1) 



2tt Jo 



sin^(x)dx 



+ / sin^(X,) 
Jvk 



Ys{Tk - r]k) 



2tt 



ds . 



and we have 



Tk 



Vk 



YsjTk - Vk) 

27r 







fTk l-Tk 


ds 










Jrik JVk 



sin^ {Xs ) drds 



< — 
- 27r 



1 fTk l-Tk 



Vk ''Vk 



2-K 



\Yr-Ys\drds. 



Finally we have the 

Lemma 6.2. Uniformly in t > 0, 

j:k^Koi;:i;:\Yr-YMrds 

T,keKo{rk - Vk) 

a.s., as M ^ oo. 

Proof: Since \Yt - Yr,J < I foi Vk < t < Tk, 
EkeKoQQ\Yr-Y,\drds 



^ / X < 2 sup (ta; - Vk) 

l^keK'^\Tk - Vk) kf^Kf 

< Svr/M 
^0, 



as M — > oo, uniformly in t. 



□ 



We are now in a position to prove the following ergodic type theorem, 
from which Theorem 13.21 will follow : 



Proposition 6.3. As t ^ oo, 

^j\in\Xs)ds^^ 

in probability. 
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Proof: We first note that 

[0,t]=B?UB^UCt, 

where 

5°= [O,t]n(u,,g^o[%,rfc]) , 
C, = [0,t]\{B^UBl). 

We have 

I f\m\Xs)ds = ^ riBo(s)sm2(X,)ds + ^ ri^i(s)sm2(X,)ds 
t Jo r Jo * t Jo ' 

+ 1 f\cAs)sm\Xs)ds. 
t Jo 

Now Ct C {s e [0,t], ly^l < M + 1}, so for each fixed M > 0, it follows 
from Lemma [53] that the last term can be made arbitrarily small, by choos- 
ing t large enough. The second term goes to zero as M — > oo, uniformly in 
t, from Lemma |6. 11 Finally the first term equals the searched limit, plus an 
error term which goes to as M ^ oo, uniformly in t, see Lemma 16.21 and 
the following fact, which follows from the combination of Lemma 16.11 and 

Lemma 15.51 : 

I 

J ^i.'Tk- Vk) 1 
in probability, as n ^ oo. □ 



We can finally proceed with the 
Proof of Theorem 13.21 Ah we have to show is that (see (|ir3l 



1 rn?/^t 1 i'2-K 1 

lim — - / sin2(X")ds = — / s\n^{x')dx' = - 



in probability. In the case v = 0, the process {(X",yj")} does not depend 
upon n, and the result follows precisely from Proposition 16.31 Now suppose 
that f 7^ 0. In that case, the result can be reformulated equivalently as 
follows. For some x € M, y 7^ 0, each t > 0, define the process {(X*, Y^), < 
s < t} as the solution of the SDE 

^ -yt X* - T 

dy,* = sin(x*) dWs , 1^0* = • 
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We need to show that 

- [\m'^{Xl)ds^— r s\n^{x')dx' 
t Jo 2tt Jo 



in probability, as t ^ oo. Note that in time t, the process starting from 
can come back near the origin. 
It is easily seen, by introducing the Markov time r|,j = inf{s > 0, jl's*] < 
M} and exploiting the strong Markov property, that 



1 

1 Jo 



t 

l{|y,*|<M}ds a.s. 



follows readily from Lemma [5.51 The rest of the argument leading to Propo- 
sition [63] is based upon limits as M ^ oo, uniformly with respect to t. It 
thus remains to check that the fact that Yq now depends upon t does not 
spoil this uniformity, which is rather obvious. □ 



Remark 6.4. This proof holds uniformly with respect to initial data 
{u,v) satisfying \v\ > a for any a > 0. 

7. Proof of Theorem 13.11 Our Theorem 13.11 now appears as a simple 
corollary of Theorem 13.21 

Proof: We first prove the Theorem for t G [0, 27r]. Then the vector P = 
is a martingale in M/^, and to prove our claim it suffices to 
show that its quadratic variation matrix converges to nt times the identity 
matrix. The diagonal elements of the matrix are 

{P('^))t = /%in2 Q(^) + COS^ Q('^))^ds = TTt 

Jo 

and we only need to compute the cross-variation 

Jo 

(7.1) = rcos(QM-g(-'))7rds. 

Jo 

Now, define (with n = vr^/^A) 
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These processes solve the stochastic differential equation 
(7.2) dUr = nVrdt, [/o" = ^° J' , 



(7.3) V,^ 



M (f') 
Po -Po 

27ri/2 



dVr = -^{sm{Ul- + Un-sm{Ut-Un)dCt 

+ J_ (cos(C/r + f/D - cos([/;" - C/f )) dSt 
(7.4) = sinC/i"dWi" 

where the process VF" is the martingale defined by Wq = and 

dWl" = TT-^^^{cos dCt - sin t/j'" d5t) . 
The quadratic variation of is 

(T^")i = vr-i /*(cos2 C/f + sin2 C/f ) vr ds = t 
Jo 

in view of the quadratic and cross variations (j3.3p of (C, 5), and this result 
does not depend on the process If^ (which follows the center of mass of 
the two particles and i^'). Thus is a standard Wiener process, and 
(C/",y") defined by (I72I)-(IZ3])-(IZ3]) satisfies the hypotheses of Theorem 
Hence, with X defined by (j4.ip . 

{pM^pi'^'))^ = /*cos(2C/;)7rds= r(l-2sin2c/;)7rds 
JO Jo 



vrt I 1 - -TT^ / sin^X^/ds' 



n2/3t7o - J 

which converges in probability to for n ^ cxd as shown in the proof of the 
theorem. 

Now we consider the process over the interval [0,47r], taking into account 
that (C, S) over the whole interval is neither a martingale nor Markov. From 
the given initial data, amplitude A and realization of (C, 5), we define a 
subsidiary set of particles M + 1 < i' < 2J\f, with initial data 

of which a.s. none coincides (modulo 27r for q) with any of the initial data 
(Qo'^^^o"'*)- Recahing that at t = 27r the law of any (for I < < Af) 
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is Gaussian with variance 21^^ (so that its probabihty density is bounded by 
l/\/27r 27r2 < l/(27r)), this ensures that, given e > 0, 

mm _ > -P\^'\<eox mm iPaV - ^27r I < e 



\<v,v'<M \<v<vi<H 

V27r 27r2 



< , < eAAV2 . 



For processes such that vainly \Pq'^ ~^2^ ^1 — ^ ^'^'^ niiuj^^j,/ \P2^ ~^2^ ^1 — 
e, the above proof holds (uniformly with respect to initial data) for the full 
set of 2J\r particles defined here. It follows that a.s. the p}'^^ —P^^^ , t G [0, 27r] , 
1 < v < 2J\f are mutually independent Wiener processes, each with variance 
7rt. This implies that pj^'^^ — Pq^\ t G [0,47r], 1 < u < Af are also mutually 
independent Wiener processes. 

For any interval [0, 2A;7r] with /c G Nq the same argument reduces the 
AA-particle problem to kM particles over [0,27r], and the proof is complete. 
□ 



Remark 7.1. Our theorem allows that Pq'^^ = Pq'^^ for some v ^ v' 
with \<v,v'<M. 

8. Perspectives. The implications o f Theorem 13.11 are twofold. First , 



for = 1, they support the observation [BeEs97l . lBeEs98al . ICEVQCll . IeIOsI ] 



that, in systems with finite M. ^ \ and j4o/ra ^ 1 (with appropriate 
scaling), the long-time behaviour of a single particle in a periodic wave field 
exhibits statistical properties approaching those of the Brownian motion. To 
complete the connection with physics literature, one must now discuss how 
the finite sums in (j2.2p approach the right hand side of (j3.2p . and how this 
implies that the solut ions of the first equation approach the solutions of the 



latter equation jSu78| . This will be discussed separately. 

Second, and this is conceptually more fundamental, for a single sample of 
the Wiener wavefield (with = oo formally), with strong spatial correla- 
tions (thanks to the single wavevector in the model), the limit A^jvcv — > oo 
leads to independence of the evolutions of all M particles, which can then be 
collectively described by the diffusion equation. While such an independence 
is often admitted without proof in physics practice, our work provides an 
explicit justification to it. We even prove a little more than usual one- or 
two-time statements, as in our limit independent in law from even 

the full evolution data {v^''\t), 0<t<T, 1 < <M - 1}. 
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This second implication is an important issue, as the acceleration of a 
passive particle is a Hamiltonian process while the diffusion equation is irre- 
versible. While the key to this irreversibility is clearly the fact that the diffu- 
sion process only relates to the momentum component of particle evolution, 
we shall further investigate the interplay of limits — > oo and A — > oo, or 
m — > 0, with stochasticity versus Hamiltonian "conservativeness" in {Q, P) 
variables in future work. 

Finally, we followed general practice in d iscussing the sto chastic acceler- 
ation problem in only one space dimension Chi79l . HaWa04 1 . This is rather 



classical, and it applies e.g. to particle motion along magnetic field lines 
in strongly magnetized plasmas ; higher-dimensional motions may call for 
different elementary models. 

APPENDIX A 

For the convenience of th e reader, we prove the following well-known 
result (see e.g. ex. IV. 3. 16 in [ReYo99f|) 



Proposition A.l. Let rj and r be two stopping times such that < 
r] < T < r] + T and Mt = jQipgdSg, where {Bt, t > 0} is a standard 
Brownian motion and {(ft, t > 0} is progressively measurable and satisfies 
Jo <^sds < kH for all t > 0. Then for all b > 0, 



sup \Mt - M^l > 6 < 2exp 



Proof: Prom the optional stopping theorem, it suffices to treat the case 
r] = 0, T = T. We have 

P( sup \Mt\ >b)<F{ sup Mt>b)+ P( inf Mt < -b) . 

0<t<T 0<t<T 0<t<T 

We estimate the first term on the right. The second one is bounded by the 
same quantity. Define for all A > 



= exp I AMt - y 



Then 



(sup Mt>b) <fI sup Mt > exp(A6 - 

0<t<T \0<t<T y 

< ex.p{X^k'^T/2 - Xb) , 
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from Doob's inequality, since {Ai^, i > 0} is a martingale with mean one. 
Optimizing the value of A, we deduce that 

from which the result follows. □ 
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